Finite-size scaling theory for explosive percolation transitions 
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The finite-size scaling (FSS) theory for continuous phase transitions has been useful in determining 
the critical behavior from the size dependent behaviors of thermodynamic quantities. When the 
phase transition is discontinuous, however, FSS approach has not been well established yet. Here, 
we develop a FSS theory for the explosive percolation transition arising in the Erdos and Renyi 
model under the Achlioptas process. A scaling function is derived based on the observed fact that 
the derivative of the curve of the order parameter at the critical point tc diverges with system size 
in a power-law manner, which is different from the conventional one based on the divergence of the 
correlation length at tc- We show that the susceptibility is also described in the same scaling form. 
Numerical simulation data for different system sizes are well collapsed on the respective scaling 
functions. 
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Explosive percolation transition (PT) occurring in a 
modified Erdos and Renyi (ER) model has attracted con- 
siderable attention in physics communities in a short pe- 
riod [H-Hl, because most studies of discontinuous phase 
transitions have been limited to equilibrium thermal sys- 
tems. The modified ER model includes an additional rule 
that discourages the formation of a giant cluster, called 
the Achlioptas process (AP) [H to the rule in the clas- 
sic ER network model [9|. According to this rule, while 
a giant cluster develops slowly, many large-size clusters 
accumulate with their numbers exceeding to those in the 
critical state. After an enforced delay, a giant cluster is 
formed at a transition point by the aggregation of the 
large-size clusters, which proceeds in an extremely short 
time. Thus, the order parameter increases suddenly. In 
analogy to the eruption of a volcano or a seismic out- 
break, this transition is called the explosive PT. 

Finite-size scaling (FSS) theory has been useful for 
characterizing phase transitions. When the phase tran- 
sition is continuous, the critical behavior of a system in 
the thermodynamic limit can be extracted from the size 
dependent behaviors of thermodynamic quantities. For 
example, the magnetization of ferromagnetic systems is 
assumed to follow a scaling form near the critical tem- 
perature Tc, 



(1) 



where L is the lattice size, f{x) is a scaling function, 
and T is temperature. (3 and v are the critical expo- 
nents associated with the magnetization and the corre- 
lation length, respectively. This scaling function was set 
up based on the fact that the correlation length, which 
diverges as T ^ \s limited to the lattice size L. When 
the phase transition is discontinuous, however, the cor- 
relation length does not diverge. Thus, a new approach 
was needed. As a successful example, a scaling function 
was derived by using the probability distribution for the 



internal energy for the g-state Potts model with g > 4 in 
two dimensions However, FSS approach for dis- 

continuous transitions arising in disordered systems has 
not been studied yet. In this Letter, we develop FSS the- 
ory for the discontinuous PT in the modified ER model 
under the AP, and obtain a scaling function, which has 
a different form from the conventional one ([1]). 

In the classic ER model, a system is composed of a 
fixed number of vertices N , which evolves as one edge is 
randomly added to it at each time step. Hereafter, time is 
defined as the number of edges added to each node. In the 
ER model under the AP, at each time step, two edges are 
randomly selected, but only one of them is actually added 
to the system to minimize the product of the sizes of the 
clusters that are connected by the potential edge. The 
ER model based on this product rule is hereafter called 
the ERPR model. While the rule in the ERPR model 
may be too artificial to fit real systems, the ERPR model 
contains an intrinsic mechanism triggering discontinuous 
phase transitions. Thus, it is meaningful to develop a 
FSS theory with this ERPR model. 

The size of the giant cluster is determined as G'Ar(t) = 
1 — sns(t), where ns{t) is the density of s-size clusters 
at time t, and the largest cluster is excluded in the sum, 
denoted by the prime in the summation. Fig[l] shows 
GN{t) versus t for different system sizes for the ERPR 
model. The curves of G]s[{t) for different TVs intersect at 
approximately one point, namely, « 0.88845(2). We 
consider the i- intercept of the tangent of the curve GAr(t) 
at tx, denoted as td{N). Then this time is calculated as 



tdiN) = tx-GN{t.)[ 
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We find that as TV 
diverges as 



oo, the derivative of Gjv(t) at 
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FIG. 1. (Color online) (a) Plot of the giant cluster size per 
node Gjv(t) versus time t for system sizes N/W^ — 128, 256, 
512, 1024 and 2048. The curves intersect at tx, which is re- 
garded as the transition point in the thermodynamic limit. 
Inset: Plot of dG'jv(t)/dt at fx versus N. The solid line is a 
reference line with a slope of 0.48. (b) The same plot with 
scaled variables G{t) = r'^G{t) and t = {t - td)r'^^^\t^, 
where r = M2{td)/M-j,{td)- Data for different system sizes are 
well collapsed in a single curve. 



with 6 ?» 0.48(1). Thus, the derivative diverges in the 
limit iV — > oo, indicating that the PT is indeed discon- 
tinuous. 

Since the dynamics of the ERFR model involves the 
selection of two edges at each time step, it may be diffi- 
cult to analytically clarify the evolution of the giant clus- 
ter. Thus, a simple model, called the cluster-aggregation 
network model, which involves the selection of one edge 
at each time step, was introduced [l2j. This cluster- 
aggregation network model is a simple representation of 
the ERFR model, because the frequency of intercluster 
connections is dominant until the percolation threshold. 
In this model, once two clusters of sizes i and j are se- 
lected with probabilities i'^ /c{t) and j"/c(i), respectively, 
where c(t) = s'^Usit), one vertex is chosen randomly 
from each selected cluster, and the two vertices are con- 
nected. When u = 1, this model reduces to the ER 
network. To put the ERPR model in this perspective, 
we measure corresponding value of w by measuring the 
probability qs (t) of an added edge being connected to a 
vertex in a cluster of size s at time t. This probability is 
given by qs{t) = s'^ns{t) / c{t) for the cluster-aggregation 
network model. Thus, by measuring qs{t) and ns(t), we 
calculate s'^ /c{t) for the ERPR model, which is denoted 
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FIG. 2. (Color online) Plot of Ks{t) versus s at various times 
for the ERPR network with system size N = 8.192 x 10^. 
The dashed line is a reference line with a slope of 0.92. The 
data correspond to t = 0.4, 0.5, 0.6, 0.7, 0.8, 0.88, U{A), and 
tz{°) (starting from the bottom). Inset: Plot of qs{td){A) 
and qs{tz){o) versus s. The solid line is a reference line with 
a slope of —1.16. 



as Ks{t). When t < td{N), Kgit) does not exhibit a 
power-law behavior, however, as t — >■ td{N), it does as 
Ks{t) - with u « 0.92 (Fig. [2). This power-law be- 
havior persists as time progresses beyond td{N). 

Since werpr ~ 0.92 is close to wer = 1, we may get a 
hint on the evolution of the ERFR model in the context 
of the analytic formulae of the ER model. The explicit 
formula for the evolution of the giant component of the 
ER model with an arbitrary initial condition was calcu- 
lated 



12l [13[. For example, the giant cluster size near 



the transition point tc 
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where M„(0) = Yl's s""n-s{0) is the initial n-th moment. 
Thus, when a system starts an evolution from ns{0) = 
5s,i, thus A/2(0) — A/3(0) ~ 1, the giant cluster size 
begins to grow from zero to a macroscopic scale continu- 
ously at the transition point tc = 1/2. 

Here, we consider a particular initial condition, in 
which M2{0) and Af3(0) depend on N. We assume that 
ris(O) follows a flat distribution, UgiO) = ng, in the range 
< s < .s,„, where s™, the size of the largest cluster at 
t ^ 0, depends on N as Sm = iV. Then, no = 2N-^'^, 
AhiO) oc iV, and M^iO) oc TV^''. Then, a FT takes place 
at tc{N) = l/2Af2(0) oc A^-'', and G{t) - r(2M2(0)i-l) 
for t > tciN) from Eq. where r = 2A'/|(0)/Af3(0) 
turns out to be in 0(1). Thus, if time t is scaled as 
t' = fAf2(0), then G{t') has a mean field behavior sim- 
ilar to the original ER case. This scaling behavior im- 
plies that t increases ~ 0{N~^) beyond tc, GN{t) does 
so by 0(1). Thus, a discontinuous phase transition oc- 
curs. Generally, if the conditions that (i) the amplitude 
A/|(0)/Af3(0) is finite and (ii) Af2(0) diverges as N ^ oo 
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are fulfilled, a discontinuous PT can occur in the ER 
model. We show below that the origin of the explosive PT 
in the ERPR model can be understood in this scheme. 

For further discussion, we reconsider the ER model 
with an initial condition ns(0) = Ss^i- As time passes, 
small-size clusters develop, and thus the cluster-size dis- 
tribution is no longer monodispcrsc. We take a cer- 
tain time t* < tc and the amplitude at that time is 
r(r) = 2Ml{t*)lNh{t*). Since the evolution of the 
ER model proceeds continuously, it holds that G^^{t) = 
r{t*)[2M2{t*){t - t*) - 1], in which the certain time t* is 
regarded as an ad hoc time origin. This relation has been 
checked numerically. We apply this logic to the evolution 
of the ERPR model below. 

We examine the amplitude M|(i)/A/3(t) as a function 
of time for the ERPR model near the transition point, 
and find that it is finite but suddenly drops to zero in the 
thermodynamic limit (Fig.[3l[a)). In the conventional ER 
model, however, Af|(t)/M3(t) is zero near the transition 
point, and thus, the transition is continuous (Fig. [Sja)). 
Interestingly, the time td{N) defined in Eq.Q locates ap- 
proximately on the verge of the sudden drop in the ampli- 
tude (Fig.|3ljb)). Thereby, we call td{N) as the triggering 
time of the explosive PT. We examine the behavior of the 
amplitude and the second moment at td{N) as a function 
of system size N more carefully. We find the amplitude 
can be fit to a power law M2 (td) / M3{td) ^ N~'^. How- 
ever, we find that k is very small and decrease slowly 
from 0.06 to 0.04, measured by successive slopes, as the 
system size increases from 10^ to 2 x 10^. Thus, one 
may regard k as zero. This suggests that the amplitude 
is independent of N in the thermodynamic limit. How- 
ever, M2(td{N)) - iV°-'^. Next, we take td{N) as the 
ad hoc time origin, and then td{N) plays a role of the 
percolation threshold as we have seen that tc{N) — >■ as 
iV -> cx) in the ER model with the flat distribution. In- 
deed, td{N) approaches t^ as in ([2]). From these results, 
we can conclude that the ERPR model satisfies the con- 
dition for a discontinuous phase transition at the ad hoc 
time origin. On the other hand, it was noticed [sj that the 
key ingredient to explosive transition was not the details 
of the dynamics during the actual explosion, but rather 
lay in the period preceding the explosion when a type of 
power keg developed. This ingredient is embodied in the 
moments of the cluster size distribution at td{N), which 
satisfies the conditions (i) and (ii). Thus, the explosive 
PT takes place. 

We derive a scaling function for the ERPR model, in- 
spired by the formula (jlj. Since o;^^^^ ^ w^^, we take 
a general form of the giant cluster size, G{t) = r{ut — 1), 
where r and u can depend on N and t. Since we take 
td{N) as the ad hoc time origin, we calculate r and u 
at td{N), and substitute t —?' t — td{N). Then, we ob- 
tain that u{td) = r^^{td)^^\t^, which is not A'hitd) as 
in the ER case. We assume that r{td) follows the same 
form Ml{td)/Ah{td) from the ER model. Thus, our ESS 
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FIG. 3. (Color online) (a) Plot of the amplitude Mi{t)/M-i{t) 
versus time t for the ERPR network (A, blue) with system 
size Afi = 2.048 x 10^ and for the ER model (□, red), (b) 
Zoom-in plot of the amplitude versus t for the ERPR network 
with different system sizes = 1.28 x 10^, 2.56 x 10*', 5.12 x 
10'' 1.024 X lO'' and 2.048 x 10^ from the top. The amplitude 
decreases more drastically as system size is increased. Various 
characteristic times are indicated on the curve for N — Ni. 
The explosion starts at td{N) and ends at tz{N). The two 
kinds of susceptibility reach their maximum at td and tc2, 
respectively, while the crossing time is independent of A'^. 



form becomes G'Ar(t) 
function and t = {t — 



= r{td)G{t), where G{x) is a scaling 
td)r-^^L. Then, 



GNit)^N-^^'fo{{t-td{N))N'^') 



(5) 



where fo{x) is a scaling function and /3/i/ = k and 
1/p = 9 + K. We confirm this scaling function with nu- 
merical data obtained for different system sizes as shown 
in Fig. [Ijb). Even though we obtain k w 0.06 — 0.04 in 
our simulation range of system size, we expect (3/9 = 
asymptotically and thus 1/D 0.48(1) and /3 = 0. We 
notice that the scaling function is in a particular form: 
Instead of following the conventional form t — ix hi its 
argument, the scaling function in ([Sj contains the form 
t — td{N), implying that the interval t — td{N) can change 
with N for a fixed t. 

We define the susceptibility in two ways. The first 
is the mean cluster size xi{t) = J2's ^^''^s{t)/ sns{t), 
which exhibits a peak at tci{N). tci{N) is larger than 
both td{N) and t^, but approaches t^ as — > cxd. The 
susceptibility at tci{N) increases with N as xi{tci{N)) ~ 



A^^ *", confirming the previous result [14|, which also be- 
haves in the same manner at td{N). Thus, Xi(t) is writ- 
ten in the form xi{t) = N^^/^ fi{{t ~ td)N'^/^), where 
fi{x) is another scaling function. Thus, ^i/D sa 0.7. 
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FIG. 4. (Color online) (a) Plot of the susceptibility xi, the 
mean cluster size, versus time in a scaled form for different 
system sizes used in Fig. [1] t = {t — fd)r~^^^^|t^, where 
r = Miitd) /M3{td)- (b) Similar plot for the susceptibility 
X2, the fluctuation of the giant cluster size, versus time in a 
scaled form. 



The second one is the fluctuation of the giant compo- 
nent sizes defined as X2{t) = Ny^{G%{t)) - {GN{t)}^- 
This quantity exhibits a peak at tc2- We find that 
X2itc2iN)) ^ N, which holds even at td(iV). Thus, 
= 1 and X2it) = Nf2{{t - td.)N^^'') with a scaling 
function /2 . The scaling behaviors are confirmed numer- 
ically in Fig. m 

To catch the physical meaning of td{N) further, we 
investigate the cluster size distribution rig (t) as a function 
of time t. In early times, ns{t) decays exponentially. As 
time passes, it exhibits a power-law behavior in small 
cluster-size region, but a hump develops in the tail region. 
We represent ns{t) as ns{t) ^ s~'^'^*'> g{s,t, N), where the 
exponent T{t) is weakly dependent on time t and ranges 
from r « 2.05(2) — 2.09(1) as time progresses beyond 
td{N). g{s,t,N) represents the shape of the hump. To 
estimate the hump size, we first plot ns{t)s'^ versus s for 
different times and define the hump size as the number 
of nodes in the shade region in FiglS] As time passes, the 
shape area increases and attains a maximum at tm{N), 
and then it reduces to zero at time t2{N). In fact, tz{N) 
is defined as the point at which the size distribution of 
finite clusters follows a power law. We obtain the value 
of T at tz as T 2.09(1) [l^. The characteristic time 
tm{N) is numerically consistent with td{N)^ indicating 
that indeed td{N) is the triggering time of the explosive 
PT, and the powder keg becomes empty at tz{N). 

In summary, we have found that the curves of the gi- 
ant component size G'Ar(i) for different system sizes Ns 
intersect at the transition point tx, and the derivative 
of GN{t) with respect to t at diverges in a power-law 
manner as N increases. This indicates that the explosive 
PT is manifestly discontinuous. We also found that the 
evolution can be regarded as kinetic cluster aggregations 
in which the connection kernel turns out to be similar to 
that for the conventional ER model near the transition 
point. Based on these facts, we have developed a finite- 
size scaling theory for for the behavior of the explosive 
PT, which has a different origin from the conventional 




FIG. 5. (Color online) Plot of s'^ns{t) versus s for different 
times t = 0.884(o), td{A), tci{o) and t^(D) for the ERPR 
model with size N = 1.024 X 10^ The hump size, the number 
of nodes contained in the shade area, becomes maximum at 
time td{N), beyond which it reduces. 



one for continuous phase transitions, and thus has a dif- 
ferent form. We have determined the critical exponents. 
The method we develop here can also be used for other 
percolating systems exhibiting discontinuous transitions, 
for example, in the interacting complex networks [16 1. 
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